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Abstract 

We have studied the structure of ^He droplets doped with magnesium atoms using density 
functional theory. We have found that the solvation properties of this system strongly depend on 
the size of the ^He droplet. For small drops, Mg resides in a deep surface state, whereas for large 
size drops it is fully solvated but radially delocalized in their interior. We have studied the 3s3p 
^Pi <— 3s^ ^So transition of the dopant, and have compared our results with experimental data 
from laser induced fluorescence (LIF). Line broadening effects due to the coupling of dynamical 
deformations of the surrounding helium with the dipole excitation of the impurity are explicitly 
taken into account. We show that the Mg radial delocalization inside large droplets may help 
reconcile the apparently contradictory solvation properties of magnesium as provided by LIF and 
electron-impact ionization experiments. The structure of ^He drops doped with two magnesium 
atoms is also studied and used to interpret the results of resonant two-photon-ionization (R2PI) 
and LIF experiments. We have found that the two solvated Mg atoms do not easily merge into a 
dimer, but rather form a weakly-bound state due to the presence of an energy barrier caused by 
the helium environment that keep them some 9.5 A apart, preventing the formation of the Mg2 
molecule. From this observation, we suggest that Mg atoms in ^He drops may form, under suitable 
conditions, a soft "foam" -like aggregate rather than coalesce into a compact metallic cluster. Our 
findings are in qualitative agreement with recent R2PI experimental evidences. We predict that, 
contrarily, Mg atoms adsorbed in '^He droplets do not form such metastable aggregates. 

PACS numbers: 36.40.-c, 32.30.Jc, 78.40.-q, 47.55.D-, 71.15.Mb, 67.40.Yv 



I. INTRODUCTION 

Optical investigations of atomic impurities in superfluid helium nanodroplets have drawn 
considerable attention in recent years,-i^ as the shifts of the electronic transition lines (atomic 
shifts) are a very useful observable to determine the location of the foreign atom attached 
to a helium drop.- In this context, the study of magnesium atoms attached to helium drops 
has unraveled an interesting and somewhat unexpected solvation behaviour as a function of 
the number (A^) of helium atoms in the drop. Diffusion Monte Carlo (DMC) calculations^ 
carried out for small drops containing a number of helium atoms up to A^ = 50 indicate 
that a Mg atom is not fully solvated in drops of sizes below A^ ~ 30. More recent quantum 
Monte Carlo calculations-'- suggest a surface Mg state for ^He clusters with up to ~ 200 
atoms. Density Functional Theory (DFT) calculations^ for ^Hcat and ^Hcat nanodroplets 
with N > 300 doped with alkaline earth atoms have shown that Mg atoms are solvated 
in their interior, in agreement with the analysis of Laser Induced Flourescence (LIF)S and 
Resonant Two-Photon-Ionization (R2PI) experiments.- LIF experiments on the absorption 
and emission spectra of Mg atoms in liquid ^He and "^He have been reported and successfully 
analyzed within a vibrating atomic bubble model, where full solvation of the impurity atom 
is assumed.— >ii A more recent experiments^ in which electron-impact ionization data from 
Mg doped ^He drops with about 10^ atoms seems to indicate that magnesium is instead 
at the surface of the droplet, in disagreement with the above mentioned LIF and R2PI 
experiments. 

There is some ambiguity associated with the notion of solvation in a helium droplet. For 
not too small droplets, one may consider that Mg is fully solvated when its position inside 
the droplet is such that its solvation energy or atomic shift do not appreciably differ from 
their asymptotic values in bulk liquid helium, as both quantities approach such limit fairly 
alongside. However, for very small drops the energy or atomic shifts of an impurity atom 
at the center of the drop may still differ appreciably from the bulk liquid values because 
there is not enough hehum to saturate these quantities. This is the case, e.g., of Mg@^He5o 
studied in Ref. |4. 

Mella et al.- have discussed how the solvation properties of magnesium are affected by the 
number of helium atoms in small ^He drops. Since DMC calculations cannot be extended 
to the very large drops involved in LIF experiments,- they could not carry out a detailed 



comparison between their calculated atomic shifts and the experiments. They also pointed 
out the sensitivity of the Mg solvation properties on apparently small differences in the 
He-Mg pair potentials available in the literature. 

The aim of the present work is to re-examine the solvation of magnesium in ^He nan- 
odroplets from the DFT perspective, extending our previous calculations^ down to drops 
with A^ ~ 50 atoms and improving the DFT approach (i) by treating the dopant as a quan- 
tal particle instead of as an external field, and (ii) by fully taking into account the coupling of 
the dipole excitation of the impurity with the dynamical deformations of the helium around 
the Mg atom. Our results confirm that full solvation of a Mg atom in ^He drops requires a 
minimum number of helium atoms, and disclose some unusual results for small drops. We 
calculate the absorption spectrum of a Mg atom attached to small and large drops, finding 
a good agreement with the experiments for the latter. We discuss in a qualitative way the 
effect of the impurity angular momentum on the electron-impact ionization yield and on 
the absorption spectrum. We also address the structure of a two-magnesium doped drop; 
the results are used to discuss the scenario proposed by Meiwes-Broer and collaborators to 
interpret their experimental results on R2PI.- 

This work is organized as follows. In Sec. 11 we briefly present our density functional 
approach, as well as results for the structure of Mg@'^HeAr drops. The method we have 
employed to obtain the atomic shifts is discussed in Sec. Ill, and applied to the case of Mg- 
doped ^He droplets in Sec. IV. In Sec. V we study how the presence of a second magnesium 
atom in the droplet may alter both the helium drop structure and the calculated atomic 
shift. A summary is presented in Sec. VI, and some technical details of our calculations are 
described in the appendices. 

II. DFT DESCRIPTION OF HELIUM NANODROPLETS 

A. Theoretical approach 

In recent years, static and time-dependent density functional methoda^^ii^ii^ have become 
increasingly popular to study inhomogeneous liquid helium systems because they provide 
an excellent compromise between accuracy and computational effort, allowing to address 
problems inaccessible to more fundamental approaches. Although DFT cannot take into 



account the atomic, discrete nature of these systems, it can nevertheless address highly 
inhomogenous helium systems at the nanoscale,— including the anisotropic deformations 
induced by atomic dopants in helium drops (see Ref. y for a recent overview on the physics 
of helium nanodroplets). 

Our starting point is the Orsay-Trento de nsity functional,— together with the Mg-He 



adiabatic ground-state potential X^T, of Ref. [l9|, here denoted as VMg-He- To check the 
sensitivity of our results to the details of different available pair potential describing the 
Mg-He interaction, we also use the sligthly less attractive potential computed in Ref. |20 . 



For the sake of comparison, we plot both potentials in Fig. [H Despite the apparently minor 
differences between these two potential curves, they cause very different solvation properties 
of Mg in small ^He drops, as we will show in the following. 

The energy of the Mg-helium system is written as a functional of the Mg wave function 



$(r) and the '^He "order parameter" \E'(r) = Jp{r), where p(r) is the ^He atomic density: 
E[^,$] = -^^- /d3r|V^(r)|2+ fdh£{p) 

2 niHe J J 

2 niMg 



+ -^/d3r|V<l>(r)p+ I Uhdh'\^v)\''VMg-He{\r-v'\)p{v'). 



In this expression, S{p) is the ^He "potential energy density".— Minimizing E under the 
constraints of a given A^ and a normalized Mg wave function, yields the ground state of 
the drop-impurity complex. To address the solvation of the Mg atom, we have found it 
convenient to minimize E subjected to the additional constraint of a fixed distance Zq 
between the centers of mass of the helium moiety and of the impurity atom which, due 
to the symmetry of the problem, can both be taken on the z axis. This is done following 
a method -borrowed from Nuclear Physics- similar to that used to describe the fission of 
rotating ^He drops.— Specifically, we minimize the expression 

E + ^[Z-Zor, (2) 

where Z is the average distance between the impurity and the center of mass of the helium 
droplet 

Z= fdr^z\<l>{r)\^-— fdr^zp{r) . (3) 

Ac is an arbitrary constant. The value Ac ~ 1000 KA^^ has been used in our calculations, 
which ensures that the desired Zq value is obtained within a 0.1 % accuracy. 



We have solved the Euler-Lagrange equations which result from the variations with re- 
spect to \E'* and $* of the constrained energy Eq. ([2]), namely 

- -^Avl/ + (^ + f/^, - X^[Z - ^o)^U = /i^ (4) 

2 niHe \op N ] 

- T> A$ + {UMg + \c{Z - Zo) z} $ = £$ , (5) 

where /x is the helium chemical potential and e is the lowest eigenvalue of the Schrodinger 
equation obeyed by the Mg atom. The above coupled equations have to be solved selfcon- 
sistently, starting from an arbitrary but reasonable choice of the unknown functions \1/ and 
$. The fields Uue and Umq are defined as 

UMg{r) = y'dVp(r') VMg^Heilr - r'l) . (6) 

In spite of the axial symmetry of the problem, we have solved the above equations in 
three-dimensional (3D) cartesian coordinates. The main reason is that these coordinates 
allow to use fast Fourier transformation techniques'^ to efficiently compute the convolution 
integrals entering the definition of S{p), i.e. the mean field helium potential and the coarse- 
grained density needed to compute the correlation term in the He density functional,— as 
well as the fields defined in Eq. IQ. 

The differential operators in Eqs. fllMSl) have been discretized using 13-point formulas 
for the derivatives. Eqs. (jMSj) have been solved employing an imaginary time method;— 



some technical details of our procedure are given in Ref. 12^. Typical calculations have been 
performed using a spatial mesh step of 0.5 A. We have checked the stability of the solutions 
against reasonable changes in the spatial mesh step. 

B. Structure and energetics of Mg-doped helium nanodroplets 

Equations (jMS]) have been solved for Ac = and several A^ values, namely A^ = 30, 
50, 100, 200, 300, 500, 1000, and 2000. This yields the ground state of the drop-impurity 
complex, and will allow us to study the atomic shift for selected cluster sizes. 

Figure [2] shows the energy of a magnesium atom in a drop, defined as the energy difference 

S^iMg) = E{Mg@''HeM) - Ei^He^) . (7) 



We compare in Fig. [2] the SN{Mg) values here obtained with those of Ref. |7|, where the 
Mg atom was treated as an infinitely massive particle -i.e., as a fixed external field acting 
on the ^He drop. The neglect of the quantum kinetic energy of the impurity overestimates 
the impurity solvation energies by quite a large amount, about 19.7 K for A^ = 50 and 18.8 
K for A^ = 2000. The value we have found for S^olMg), —18.4 K, compares well with the 
DMC result of Ref. lO (—21 K), showing that DFT performs quite well for small clusters,- 
far from the regime for which it was parametrized. The DMC energy found for the same 
system is ~ —168.2 K,^ whereas our DFT result is ~ —157.0 K, and the "asymptotic" DMC 
value for S2oo{Mg), —33.1 K, compares well with the DFT value, still far from the limit 
value corresponding to a very large helium drop (see Fig. [2]). 

The solvation properties of the Mg atom are determined by a delicate balance between the 
different energy terms -surface, bulk and helium- impurity- in Eq. ([1]), whose contribution 
is hard to disentangle and depends on the number of atoms in the drop, as shown by DMC 
and DFT calculations. To gain more insight into the solvation process of magnesium in 
small ^He droplets, we have computed the energy of the doped droplet as a function of the 
impurity position Zq. 

The bottom panel in Fig. [3]shows E{Zq) for N = 50, as computed using the two different 
Mg-He pair interactions shown in Fig. [1] (the discussion of the top panel is postponed until 
Sec. IV. A). It can be seen that E{Zq) in both cases displays two local minima. In one case 
(i.e. for the sligtly less attractive potential in Fig. [1]) a "surface" state for the Mg atom 
is energetically preferred, while in the other case the impurity prefers to sit in the interior 
to the droplet (although not exactly at its center). These results are in agreement with 
the DMC calculations of Ref. 14, where it has been shown that a bimodal distributions for 
the Mg radial probability density function with respect to the center-of-mass of the helium 
moiety appears for A^ < 30.- More recent DMC calculations carried out up to A^ ~ 200 
drops^ seem to point out that Mg is always in a surface state, although somewhat beneath 
the drop surface. Our DFT calculations yield that Mg is already solvated for A^ = 200. 

For both pair potentials, the bottom panel of Fig. [3] shows that the two local minima 
in E{Zq) are separated by an energy barrier of about 1 K height, allowing the impurity to 
temporarily visit, even at the experimental temperature T ~ 0.4 K, the less energetically 
favored site. This causes changes in the total energy of the system by less that 1 %, but has 
a large effect on the value of the atomic shift. We will address this important issue, and its 



consequencies on the computed spectral properties of Mg@HeAr in the next Sections. 

FigureHlshows the hehum configurations for the four stationary points displayed in Fig. [3l 
namely those corresponding to Zq = 0, 2, 4, and 6 A. Although the preference for the surface 
or the solvated state depends on the He-Mg pair potential used in the calculations, similarly 
to the case of other alkali atoms/ this does not seem to be the case for the stationary points 
at ^0 = 2 and 4 A that are present in both curves. We have compared in Fig. [5] the density 
profiles along the z axis for the Zq = 2(4) A configuration with the profile of the pure ^He 
drop, finding that the appearence of these stationary points is related to the position of the 
density peak in the first helium solvation shell with respect to a maximum(minimum) of the 
density of the pure drop. We see that a minimum(maximum) in the energy is associated 
with a constructive(destructive) interference in the oscillation pattern of the He density. We 
are lead to conclude that the interplay between the density oscillations already present in 
pure drops, and the solvation shells generated by the impurity plays an important role in 
the solvation properties of the Mg atom. This effect is also present, although to a lesser 
extent, in Ca-doped ^He nanodroplets.— 

Eventually, for larger drops Mg becomes fully solvated. We have found that this is the 
case whichever of these two potentials we use (see for instance the bottom panel of Fig. [6]). 
For this reason, the results we discuss in the following have been obtained with the Mg-He 



pair potential of Ref. [191, unless differently stated. 

When the Mg atom is fully solvated, e.g. for A^ = 1000, we have found that E{Zq) grows 
monotonously as Zq increases (i.e. as the Mg atom approaches the droplet surface). This 
is shown in the bottom panel of Fig. O To better understand how E{Zq) depends on A^, 
we have plotted in the top panel of Fig. [6] the energy of the A^ = 1000 and 2000 doped 
drops, referred to their equilibrium value, as a function of the distance from the dividing 
surface, i.e. the radius R1/2 at which the density of the pure drop equals pfe/2, pb being the 
liquid density value (-R1/2 = tqN^^^, with tq = 2.22 A). These radii are 22.2 and 28.0 A for 
A^ = 1000 and 2000, respectively. 

The top panel of Fig. [H] shows that most of the change in AE takes place in the 15 
A outer region of the drop, irrespective of its size. As a consequence of the flatness of 
E{Zo), the magnesium atom is very delocalized in the radial direction even in relatively 
small drops. This delocalization might affect the absorption spectrum of the attached Mg 
atom, and also must be explicitly considered in the interpretation of the electron-impact 
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yield experiments.— We address this issue in Sec. IV. 

III. EXCITATION SPECTRUM OF A MG ATOM ATTACHED TO A "^HE DROP 

Lax methods^ offers a realistic way to study the absorption spectrum of a foreign atom 
embedded in liquid drops. It makes use of the Franck-Condon principle within a semiclassical 
approach, and it has been widely employed to study the absorption spectrum of atomic 
dopants attached to ''He drops, see e.g. Ref. l25| and references therein. The method is 



usually applied in conjunction with the diatomics- in- molecules theory,— which means that 
the atom-drop complex is treated as a diatomic molecule, where the helium moiety plays a 
role of the other atom. 

In the original formalism, to obtain the line shape one has to carry out an average on the 
possible initial states of the system that may be thermally populated. Usually, this average 
is not needed for helium drops, as their temperature, about 0.4 K,— is much smaller than the 
vibrational excitation energies of the Mg atom in the mean field represented by the second 
of Eqs. ([6]). In small helium droplets thermal effects can show up in the Mg absorption 
spectrum due to the high mobility of the atom. For large drops, thermal motion plays a 
minor role, as the Mg atom hardly gets close enough to the drop surface to have some effect 
on the line shape. In this case, however, dynamical deformations of the cavity around the 
impurity may be relevant.—*^ 

The line shape for electronic transitions from the ground state (gs) to the excited state 
(ex) in a condensed phase system can be written as 

I{uj) (X Jdt e-'"*(^3^|Dj, e^^- Dg, e-ir««»|^5-) , (8) 

where Dg^ is the matrix element of the electric dipole operator, |\I^^'^) is the ground state 
of the system, and Tigs and Tigz are the Hamiltonians that describe the ground and excited 
states of the system respectively. 

I{uj) is evaluated using the Born-Oppenheimer approximation and the Franck-Condon 
principle, whereby the heavy nuclei do not change their positions or momenta during the 
electronic transition. If the excited electron belongs to the impurity, the helium cluster 
remains frozen, so that the relevant coordinate is the relative position r between the cluster 
and the impurity. This principle amounts to assuming that Dg^ is independent of the 



nuclear coordinates. Taking into account that e n'^3'>\'$9s^ _ ^ itujgs^^gs^ ^^^ projecting on 
eigenstates of the orbital angular momentum of the excited electron \m) one obtains 



/M oc \D. 



|2 



gel 



^ / dt e-"^^+^x)i I d"[a] I d\ ^f (r, [a])* e^^™ ('''["l) ^f (r, [a]) , (9) 



where huj^ and 'ip^i^i W\) ^'^^ the energy and the wave function of the ro- vibrational ground 
state of the frozen helium-impurity system, and H^{r, [a]) is the ro- vibrational excited 
Hamiltonian with potential energy ^^(r, [a]) determined by the electronic energy eigen- 
value, as obtained for a p <— s transition. At this point, we have introduced the variables [a] 
to represent the degrees of freedom needed to describe possible deformations of the system, 
corresponding to the zero point oscillations of the helium bubble around the impurity. If 
this effect is neglected, the deformation parameters [a] are dropped and the ground state 
wave function tpx coincides with the Mg wavefunction $ found by solving Eq. ([5]). 

If the relevant excited states for the transition have large quantum numbers, they can 
be treated as approximately classical using the averaged energy hu^ ^ V-^i'^, [c^]) which is 
independent of z/. In this case we obtain the expression 

JM ex E/d"H/d^r|V^f(r,[a])r<5(c. + c.f-Vr(r,N)A) 



hfd-\a]f d^r ^^^MM^llL no) 



where flm{^) is the surface defined by the equation u + uj^ — V^{r, [a])/h = 0. If the atom 
is in bulk liquid helium, or at the center of the drop, the problem has spherical symmetry 
and the above equation reduces to 



47fhJ2 /d"[a] 



[r^f(r,H)]^ 



dV^--{r, [a])/dr 



(11) 



where rm{uj) is the root of the equation u + Ux — Krf (^) l*^])/^ = 0- ^^ ^^e non-spherical 
case, we have evaluated /(cu) using the first expression in Eq. (TTOl) . 

The potential energy surfaces V^{r, [a]) needed to carry out the calculation of the atomic 
shifts have been obtained in the pairwise sum approximation, using the Vuir) and Vj]{r) Mg- 
He adiabatic potentials from Ref. \4- In cartesian coordinates, and assuming that the He- 
impurity spin-orbit interaction is negligible for magnesium, the eigenvalues of the symmetric 
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matrix 

t/.,(r, [a]) =Jdh'p{r' + r, [a]) |vn(r')5., + [V^{r') - Vn{r')]^^ (12) 

are the V^{r, [a]) potentials which define the potential energy surfaces (PES) as a function 
of the distance between the center-of-mass of the droplet and that of the impurity. 

For spherical geometries, Eq. (fT2l) is diagonal with matrix elements (in spherical coordi- 
nates) 



Xi{r, [a]) = U^i{r, [a]) = 2n f f r'^ sin 9'd9'dr'p{Vr'^ + r^ + 2r'r cos 9', [a]) 



X Vu{r') + [V^{r') - Vuir')] 



^{Sii + Si2) sin^ 9' + 6i3 cos^ 9' 



(13) 



In this case, two of the PES are degenerate, namely Ai(r, [a]) = A2(r, [a]) ^ Xs^r, [a]).— 
This holds true for r 7^ 0, and it is relevant when we take into account the delocalization of 
the impurity inside the bubble due to its quantum motion. Otherwise, since at r = all the 
Aj coincide, they are threefold degenerate. 

IV. RESULTS FOR THE ABSORPTION SPECTRUM OF MAGNESIUM ATOMS 

The problem of obtaining the atomic shift of magnesium in a helium drop has been 
thus reduced to that of the dopant in the 3D trapping potentials corresponding to the 
ground state and P excited states. We consider first the [a] = case (i.e. no zero-point 
deformations of the helium cavity hosting the impurity). The general situation, in particular 
the homogeneous width calculation, is presented later on in this Section. 

If [a] = 0, the model is expected to yield at most the energies of the atomic transitions, 
but not the line shapes since the impurity-droplet excitation interactions as well as inho- 
mogeneous broadening resulting from droplet size distributions, laser line width and similar 
effects are not included. These limitations are often overcome by introducing line shape 
functions or convoluting the calculated lines with some effective line profiles.— i^^ We discuss 
here two illustrative examples, namely the atomic shift of magnesium in A^ = 50 and 1000 
nanodroplets. The homogeneous width is calculated in Subsection B for the A^ = 1000 
droplet. 

For Mg@^He5o, the calculated shifts at Zq = (spherical configuration), 2 and 6 A 
are 500, 450 and 281 cm~^, respectively (281.3, 281.7 and 283.0 nm wavelengths). No 
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experimental information is available for such a small drop. Contrarily, the absorption 
spectrum of the 3s3p ^Pi ^- 3s^ ^Sq transition of Mg atoms attached to large helium drops 
has been measured,- displaying a broad peak strongly blue-shifted from its position in the 
gas-phase. This spectrum is remarkably close to the one obtained by Moriwaki and Morita^- 
in bulk liquid helium, and hence it has been concluded that Mg is in the interior of the ^He 
droplet. We want to point out that, while the absorption line in liquid helium was attributed 
to a single broad peak of energy 281.5 nm (35 524 cm~^), i.e., a shift of 474 cm~^,^- in large 
drops a similar line profile^ was fitted by two Gaussians centered at 35 358 and 35 693 cm~^ 
(282.8 and 280.2 nm wavelengths) respectively, i.e., shifted 307 cm~^ and 642 cm~^ from the 
gas-phase line.- The origin of the two peaks was attributed to the splitting of the degenerate 
n state by dynamical quadrupole deformations of the cavity surrounding the dopant, since 
this argument had qualitatively explained similar doubly-shaped D2 excitation spectra of Rb 
and Cs atoms in liquid "^He due to a quadrupole oscillation of the hehum bubble (dynamic 
Jahn- Teller effect).— 

LIF experiments on the heavier alkaline earth Ca and Sr in large ^He droplets^ have 
disclosed the existence of strong blue-shifted, broad peaks with no apparent structure, al- 
though it cannot be discarded that this broad line could be a superposition of unresolved 
peaks. The same happens for Ba.— The surface location of Ca, Sr, and Ba in these drops 
has been further confirmed by DFT calculations.-i2^ It is also interesting to recall that LIF 
experiments on Ca atoms in liquid *^He and ^He have found a broad line in the region of the 
4s4p ^Pi <— 4s^ ^So transition with no apparent splitting,— contrarily to the case of Mg. 

Since Mg is fully solvated in the A^ = 1000 drop, the calculated atomic shift Au may be 
sensibly compared with the experimental data where drops with A^ in the 10'^ — 10^ range 
are studied. We have obtained Au = 659.0 cm~^ (280.0 nm wavelength); this peak nearly 
corresponds to the Gaussian that takes most of the intensity of the absorption line (about 
87 %).- We have carried out a detailed analysis for this drop, determining the equilibrium 
structure of MgQ'^Heiooo as a function of Zq, and have used it to evaluate Au. The results 
are displayed in Table [H showing the actual sensitivity of the absorption spectrum to the 
Mg atom environment. 

The impurity-drop excitations will determine the homogeneous width of the spectral line, 
and the population of excited states may be relevant given the limit temperature attained 
by the droplets.-*^"^ In this context, the relevant excitation modes of the helium bubble 
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are radial oscillations of monopole type (breathing modes), and multipole shape oscillations 
about the equilibrium configuration, as well as displacements of the helium bubble inside 
the droplet. We will address these issues in the following Subsections. 

A. Thermal motion and angular momentum effects 

To describe the displacement of the helium bubble inside the droplet, we have fitted 
the E{Zq) curve of the Mg@^Heiooo system to a parabola, and have obtained the excitation 
energy huj for this 3D isotropic harmonic oscillator. The hydrodynamic mass of the impurity 
atom has been estimated by its bulk liquid helium value, M* ~ 40 a.u., obtained by the 
method outlined in Appendix A, Eq. (]A23p . We find hu = 0.1 K, indicating that thermal 
motion, i.e., the population of the excited states of the "mean field" E{Zq) is important 
at the experimental temperature T = 0.4 K, and may produce observable effects in the 
absorption spectrum and the electron-impact ionization yield. 

To describe in more detail the delocalization of the Mg atom inside the drop, we have 
used an effective Hamiltonian where we interpret Zq as the radial distance R between the 
impurity and center of mass of the helium moiety, and E{Zq) as the "potential energy" 
V{R) associated with this new degree of freedom of the impurity in the drop. Namely, 

n = —- + V(R) = -^ + , + V(R) , (14) 

2M* ^ ' 2M* 2M*R^ ^ v ; , K ) 

where M* is the Mg hydrodynamic mass. In the canonical ensemble, the total probability 
distribution W a.s a. function of R can be written as 



W[R) =Q^'J di?'i?'2|dfi'5:(V^„,^|R')e-^(^')/^-^(RVn£m) 
= Q-' ri?'Mi?'^(2£+ l)e-^-/^^^|0„,(i?')l' , 



where the partition function is defined as Q = Tr ( e '^/^bT\ _ ^^^(2£ + l)e ^rii/ksT ^ being 
kB the Boltzmann constant. The radial probability density w is 

w{R) = ^ = Q-'R'Y.i'^l+l)e-^-^/'-''\<^rAR)? ■ (16) 

This expression has been evaluated for A^ = 50 and 1000 by solving the Schrodinger equation 
for the Hamiltonian Eq. ^\M to obtain the orbitals 0„£ and eigenenergies En^. For larger 
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N values, we have used the semiclassical approximation Eni — > 2m* ~^ Ves{R)-, where the 
effective potential is 

Integrating pr in phase space, we obtain for the probability density 

1 h^^{i + iy 



w{R) = Q-'R'e^v (-^) E(2^+ l)e^P 



keT 2M*i?2 



(18) 



1 h^i(i+i) 



with the normalization Q = J^ dRR"^ exp [ — -jjjrj J^tC^^ + 1) 6xp ^ ^ 2M*ij2 

Lacking a better choice, we have weighted any possible angular momentum value with a 
Boltzmann energy factor. It has been shown^S that some of the angular momentum deposited 
in the droplet during the pickup process may be kept in the impurity atom, resulting in a 
different angular momentum distribution than the Boltzmann one. This could yield that 
some Mg atoms are actually closer to the drop surface. 

If the angular momentum associated with the motion of the magnesium atom -whose 
"radius" is ~ 5 A, see Fig. [1]- is such that Mg can be some 10 A beneath the drop surface, 
the shift of the absorption line would be hardly distinguishable from that of the totally 
solvated case -as seen in Table [H At the same time, the electron-energy dependence of the 
Mg"*" yield observed in electron impact ionization experiments^^ (and which was considered 
as an evidence of a surface location of Mg atoms on ^He droplets) could indeed be due to 
Penning ionization of the impurity in a collision with a metastable He* atom that occupies 
a surface bubble state in the drop, instead of being due to the transfer of a positive hole 
(He"'") to the Mg atom, which is the primary ionization mechanism when the impurity is 
very attractive and resides in the deep bulk of the droplet. 

We show in the top panel of Fig. [3] the probability densities w{R) at T = 0.4 K corre- 
sponding to the configurations displayed in the bottom panel.^ Similarly, the top panel of 
Fig. [6] shows that for Mg@^Heiooo and Mg@^He2ooo; if thermal motion is taken into account 
and the impurity retains some of the pick-up angular momentum, the maximum density 
probability of Mg is at ~ 15 A beneath the drop surface in both cases. To obtain it, we 
have taken for M* the bulk value 40 a.u. As seen from Table [U the absorption line shift 
changes by a small 2% with respect to the R = configuration. The values of the angular 
momentum corresponding to these maximum density probabilities are (L^)^/^ ~ 9^ and 
~ 10^, respectively.*^ For a A^ = 10 000 drop,— whose radius is R1/2 = 47.8 A, we have 
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extrapolated inwards the AE curves of the calculated A^ = 1000 and 2000 drops, and have 
obtained from it the probability distribution displayed in Fig. [61 Its maximum is at ~ 18.5 
A beneath the surface, with (L^)-*^/^ ~ 28^. 

Finally, using the effective potential of Eq. flTT]) we have determined that for A^ = 
1000(2000), a Mg atom with d. ~ 2Ah{AAfi) is in an "equilibrium position" some 10 A 
beneath the drop surface. For the A^ = 10 000 drop, this value is £ ~ 148/i. These values 
look reasonable, and Mg atoms holding this angular momentum or larger might thus be 
the origin of the primary electron-collision ionization yield by the Penning process, without 
questioning the conclusion drawn from LIF experiments that magnesium is fully solvated in 
"^He drops. 

B. Homogeneous width from shape deformations of the helium bubble 

We have shown that for large drops, the Mg atom is fully solvated and its thermal motion 
only produces small changes in the absorption shift. This allows us to decouple the effect of 
the translational motion of the helium bubble on the absorption line, from that of its shape 
fluctuation. Moreover, we can address shape fluctuations in the much simpler spherically 
symetric ground state, when magnesium is located at the center of the drop. 

To quantify the effect of these fluctuations, we have first used the spherical cap model^ to 
estimate the excitation energies of the helium bubble around the impurity in liquid helium. 
To this end, we have fitted the VMg-He potential to a Lennard- Jones potential with depth 
e = 7 K at a minimum distance Vmin = 5 A. Minimizing the total energy within this 
model yields a configuration with an equihbrium radius of Rq = 0.97 x 2^^/^rmm, that we 
approximate by Rq = 2^^l'°rmin to obtain the excitation energies. 

Deformations of the ^He around the Mg atom are modeled as^*^ 



R{Sl) = Ra 



oo A 
A=2/^=-A 



(19) 



where R is the radius of the bubble cavity hosting the solvated Mg atom, Vt represents the 
solid angle variables {9, (p), Zxf^{Q) is a real spherical harmonic, and ax^ is the amplitude of 
the A/i— multipole deformation. The condition R^ = Ro{l — -^J^x/i Ic^A/^P) ensures the con- 
servation of the number of particles up to second order in axfj,- The dipole mode amplitude 
aifj, is absent since, for an incompressible fluid, it corresponds to a translation of the bubble, 
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and this has been considered in the previous Subsection. 

If ^[/^(fi)] represents the bubble surface and a the surface tension, the energy for a large 
drop can be written as 

12 



/■oo 

E = aS[R(n)]+pbe / d^r 

jR(n) 



R{n) [\ r 
1„ , ,o 1 °^ ^ 



-2 



^ E,, + -Co\ao\' + 2 E Ca E I^AmI' , (20) 

A=2 /i=— A 

where Cq = 87rcri?Q (1 + 12Aa) and Ca = cr-Roll-^ ~ l)!-^ + 2) + 6A^] are the stiffness pa- 
rameters, and A^ = pbs2~^^^rminO'~^ is the impurity- He solvation parameter.— The mass 
parameters are Bq = 4:7ipi,mHeRo and Bx = pbrriHeRo/i^ + l),"^^ and the excitation energies 



are determined from Hlux = hJCx/Bx, yielding hujQ = 10.2 K for the breathing mode, and 
hu!2 = 9.8 K for the quadrupole mode. Given the droplet temperature of 0.4 K, we conclude 
that only the ground state is populated. The mean amplitude of the shape oscillations is 
estimated from the variance 7a = h^ {BxCx)~ /2, giving 70 = 0.03 and 72 = 0.15. This 
model thus yields that the bubble can experience monopole oscillations of ~3% amplitude, 



and quadrupole deformations of 72A/3/47r ~ 8% amplitude. Amplitudes of this order have 
been determined within the atomic bubble model for Cesium atoms in liquid helium.— Since 
their effect in the absorption spectrum is expected to be relevant, we have undertanken a 
more refined calculation within DFT taking Mg@^Heiooo as a case study. 

For_helium droplets, we have described bubble deformations in a way similar as in Refs. 
, namely, if po{r) is the helium spherical ground state density, deformations are 
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introduced as p(r,t) = pQ[R{r,t)] /C ^, with 

00 A 

R{r, t) = r + aoit) + E E «Am W^A^(r) , (21) 

A=2/^=-A 

where the real spherical harmonics are normalized as {Zx^lZy^') = 2i+T^AA"^/i^' for conve- 
nience, and the normalization /C = N~^ J d^rpQ[R{r,t)] ensures particle number conserva- 
tion. If the Mg wave function follows adiabatically the helium density deformation, it can 
be shown that to second order in axfi, the total energy of the system can be written as 

E{[a],[a]) ~ Egs + lM*al + 27TEi,'^al 

+ i:EJiMJ"5„.^4='„5„}, (22) 
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being Egg the ground state energy, M^ the hydrodynamic mass asociated with the A mode, 

(2) 

and E\ the second derivative of the total energy with respect to axf^. This equation 
represents the Hamiltonian of a set of uncoupled harmonic oscillators, whose quantization 
yields a ground state to whose energy each mode contributes with e^^ = I^^a? with ux = 
, and a ground state wave function 



(2A+1)M* 



na]) = [-;^) e - °n(;^ e .« ^, (23) 

where a| = Y.^=~.x Q^a^- Details are given in Appendix A. 

We have computed the hydrodynamic masses assuming that the drop is large enough 
to use Eq. (lA23p . We have obtained Mq = Ib.OniHe + 0.28mMg ~ 66.7 a.u. and M2 = 
1.9m//e+0.56mA/g ~ 21.0 a.u. In actual calculations, instead of using Eq. flAlOp . the energies 



(2) (2) 

Eq and E2 have been numerically obtained by computing the total energy of the system 

for different small values of ao and 02- This has been carried out by numerically introducing 

the desired deformation parameter into the ground state density and renormalizating it, 

solving next the Schrodinger equation for the Mg atom [Eq. (JS])] to determine the ground 

state of the impurity, and computing the total energy of the system from Eq. ([1]). Fitting 

these curves to a parabola, we have obtained Eq = 49.7 KA~^ and E2 = 16.8 KA~^. We 

have then calculated the ground state energies huJx/2 and deformation mean amplitudes 7^, 

obtaining hujQ/2 = 10.6 K and huj2/2 = 6.3 K, with mean amplitudes 70 = 0.18 A (~3.7%) 

and 72 = 0.42 A (-8.5%). 

To quantitatively determine the effect of these deformations on the absorption spectrum, 
we have developed Eq. (TTO!) to first order in the deformation parameters, and have explicitly 
shown that to this order, only the breathing and quadrupole modes affect the dipole ab- 
sorption spectrum. The details are given in Appendix B, where we show that the breathing 
mode affects the shift and shape of the line, whereas quadrupole modes only affect the shift. 

Consequently, we restrict in Eq. ( ITT]) the deformation parameters needed to properly 
describe the homogeneous broadening of the absorption dipole line, namely d[a] — ^ dao d^a2 
and ^x(^5 H) ^ V^(q^O) C(2)^{r, ao), being $(r, ao) the wave function of the Mg atom for a 
given ao value, and compute the spectrum as 



I{uj) oc 47r /daod^a2|V^(ao,a2)P5I / ^^ 1^ *('^'"o)P 



x6 



UJ + tuf (ao) - J^^m{f^ tto, ["2]) 
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47r / daod^a2\ip{aQ,a2)\'^h^ 



[r $(r,ao)]^ 



(24) 

r=rm{u!) 



dV^{r,ao, [02]) /dr 

where [^2] = {"2-2, "2-1, "20, «2i, "22}, V^{r,ao, [02]) are the eigenvalues Am(r, [a]) of the 
excited potential matrix of Eq. flC4p . and rm{uj) is the root of the equation u + ujx{o:o) — 

Expression ([23]) has been integrated using a Monte Carlo method. We have sorted 
M = 10^ sets of values [a]* = {oq, a2_2, «2-i; '^20' '^21; '^22} using the square of the wave 
function of Eq. (!23|) as probability density. Next, for each set we have found the eigenvalues 
V^{r, ag, [02]*) of the Uij matrix that define the potential energy surfaces and have used a 
trapezoidal rule to evaluate the integrals 



Uu^, [af) =471 I dr \r <l>{r,al)\'6 cu + cuf K) - V-(r,a^, H^ 



(25) 



h 

using a discretized representation of the delta function.-^ Finally, we have obtained the 
spectrum as 

1 M 

^H«T7EE^-(^J«n- (26) 

Figure [3 shows the absorption spectrum of one Mg atom attached to ^Heiooo in the 
vicinity of the 3s3p ^Pi ^- 3s^ ^Sq transition when homegeneous broadening is considered. 
We have decomposed the absorption line into its three components, the higher frequency 
component being the S one. The starred vertical line represents the gas-phase transition, 
and the experimental curve, adapted from Ref. ISi, has been vertically offset for clarity. Also 
shown is the absorption spectrum obtained by neglecting homogeneous broadening (hatched 
region). This figure shows that the both the energy and width of the absorption peak that 
takes most of the experimental intensity are correctly described by our calculations. 

V. TWO MAGNESIUM ATOMS ATTACHED TO A ^HE DROP 

The attachment of magnesium atoms in ^He droplets has been recently addressed using 
resonant two-photon-ionization.-^ In particular, the authors of Ref. M have obtained the 
absorption spectrum for drops doped with different, selected numbers of Mg atoms. From 
their measurements it appears that two main features contribute to the observed line shapes, 
one peaked at about 279 nm, and another at about 282 nm. This is in agreement with the 



results of Refs. 
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(we recall that actually, the two peaks were not resolved by the authors 
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of the bulk liquid experiment). The structure at 282 nm, however, only appears if the 
droplet contains more than one Mg atom. Thus the two-peak structure cannot be due to 
the splitting of the absorption line due to dynamical quadrupole deformations of the helium 
bubble around the impurity, as previously believed. We have indeed shown in the previous 
Section that this coupling only produces a broad peak, in good agreement with the results 
of Ref. y for helium drops containing just one Mg atom. 

Another interesting observation reported in Ref. |9| is that their experimental results for 
multi-atom doped ^He droplets are not consistent with the formation of compact, metallic 
Mg clusters inside the ^He droplet. The magnesium atoms in the droplet appear instead to 
be relatively isolated from each other, showing only a weak interaction and leading to the 
282 nm shift in the observation. 

To confirm this scenario and find an explanation for the origin of the low energy com- 
ponent in the absorption peak, we have carried out DFT calculations to determine the 
structure of a two-magnesium doped ^He drop. Our goal is to verify whether the helium 
density oscillation around a magnesium atom may result in an energy barrier preventing 
the Mg atoms from merging into a Mg2 dimer, as suggested by Przystawik et al.^ To obtain 
the structure of two Mg atom in a ^He drop, we have minimized the energy of the system 
written as 

2 Tflij^ J J 

+ -^/d3r{|V<l>i(r)|^ + |V<l>2(r)p} 

2mMgJ 

+ ||d3rdV{|$i(r)p + |$2(r)p}l^M,_^e(|r-r'|)p(r') 

+ / /d3rdV|<l>i(r)pV^M,-Af,(|r-r'|)|<l>2(r')|' , (27) 
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and the other ingredients 



where Vug-Mgi]'^ ~ ^"'1) is the Mg-Mg pair potential of Ref. 
have the same meaning as in Eq. ([1]). 

There are at least two additional effects which are not considered when modeling the Mg- 
Mg interaction via the pair-potential in vacuum, as implied in the above expression. The first 
is due to three-body (and higher) correlation effects involving the ^He atoms surrounding 
the Mg pair: these should exert an additional, albeit small, screening effect due to He 
polarization, which is expected to reduce the absolute value of the dispersion coefficients in 
the long-range part of the Mg-Mg interaction. The second is a possible reduction of the Mg 

19 



atom polarizability due to the presence of the surrounding ^He cavity which, because of the 
repulsive character of the electron-He interaction, should make the electronic distribution 
of the impurity atom slightly "stiffer", thus reducing further the values of the dispersion 
coefficients in the Mg-Mg pair interactions. Although in principle these effects might reduce 
the net interaction between a pair of Mg atoms embedded in liquid ^He, in practice in the 
present system they are indeed very small. The correction to the leading term of the long- 
range dispersion interaction, —C^/r^, due to three-body correlation effects can be written 
to first order— as —€^{1 — 27ma/3)/r^, a being the static polarizability of the host fiuid 
{aHe = 1.39 Og). Such correction is of the order of only 1% in our case. To estimate the 
change in the Mg atomic polarizability due to the surrounding He, we computed, using ab- 
initio pseudopotential calculations, the (static) polarizability of a Mg atom in the presence 
of an effective (mainly repulsive) potential acting on the Mg valence electrons due to the 
presence of the surrounding He. The effective interaction is derived from the equilibrium 
shape of the "^He bubble hosting the Mg atom, as predicted by our DFT calculations, and 
assuming a (local) electron-He density-dependent interaction which was proposed by Cole 
et al.— We find a very small change in the static atomic polarizability a of Mg. Assuming 
that, roughly, Cq oc a^ we find a reduction of the Cq coefficient of about 1-2%. 

The minimisation of the total energy functional written above under the constraint of 
a given number of helium atoms and normalized Mg ground state wave functions should 
in principle yield the equilibrium configuration of the system. In practice, depending on 
the initial configuration, we have found several local minima, whose origin is again the 
"interference" of the He solvation shells around the Mg atoms. We have found three such 
metastable configurations for (Mg+Mg)@^Heiooo if "we start the minimization procedure 
with one Mg atom near the center of the droplet, and the other placed off center, at some 
distance from the first. They are displayed in Fig. [HI The energy difference between the 
innermost (Mg-Mg distance d = 9.3 A) and the outermost {d = 18.5 A) configurations is 
12.5 K. The energy of the d = 9.3 A configuration is sensibly that of the configuration 
specularly symmetric about the z = plane (-5581.4 K) also shown in the figure. It is worth 
noticing that, since R1/2 = 22.2 A and the "radius" of Mg is ~ 5 A, only the upper left 
corner configuration has the Mg impurity in a surface state. 

Figure [9] shows two density profiles of the symmetric configuration obtained along the 
2;-axis (solid line) and the x- or y-axis (dashed line). It shows a relatively high density helium 
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ring around the two-bubble waist, clearly visible in Fig. [HI where the local density is almost 
three times the bulk liquid "^He density, and that prevents the collapse of the two-bubble 
configuration. The Mg wave functions are peaked at ~ ±4.75 A, and very narrow. This 
justifies a posteriori the assumptions we have made to write the total energy of the system, 
Eq. m- 

Our results confirm the existence of the energy barrier suggested in Ref. |9|, that prevents 
the two Mg atoms from coming closer than some 9 A, and thus hindering, at least temporar- 
ily, the formation of the Mg2 dimer. This barrier is shown in Fig. [10] as a function of the 
Mg-Mg distance d, which is kept fixed in a constrained energy minimization. Note that the 
energy of the Mg+Mg system increases as d does because the two Mg atoms in a drop form 
a state more bound than that of a pure drop with the impurities well apart. 

Since the height of the barrier is larger than the experimental temperatures, T ~ OAK, 
two solvated Mg atoms will not easily merge into a dimer, but rather form of a metastable 
weakly bound state. Based on these finding, we suggest that several Mg atoms solvated inside 
^He drops might form a sparse, weakly-bound "foam" -like aggregate rather than coalesce 
into a more tightly bound metallic cluster. Partial coagulation of impurities was already 
invoked by Toennies and coworkers to explain their experimental findings for the successive 
capture of forcing atoms and molecules in helium clusters^ (see also Ref. 
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for the case of 



bulk liquid helium). Very recently, a kind of "quantum gel" has been predicted to be formed 
in ^He drops doped with neon atoms.— Although some degree of mutual isolation between 
foreign atoms is expected in the case of strongly attractive impurities (like those studied in 
the two cases mentioned above), where they are kept apart by the presence a solid ^He layer 
coating the impurity,^^ our calculations show that this effect is possible even for relatively 
weakly attractive impurities like Mg, where such solid-like ^He layer is absent. 
One may estimate the mean life of the metastable state as 



27r. 



^^'' exp[AU/{kBT)], (28) 



\U"id,,) 

where fi\^g = Mljg/2 ~ 20 a.u. is the hydrodynamic reduced mass of the Mg-fMg system 
and AU is the barrier height. From Fig. [lOlwe have that U"{deq) ~ 40 K A^^. This yields a 
mean life of a few nanoseconds, which is about five to six orders of magnitude smaller than 
the time needed for its experimental detection.-'-^ ^]^q mean life becomes increasingly large 
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as the relative angular momentum L deposited into the two Mg system increases. Writing 

one obtains the L— dependent energy barriers displayed in Fig. [TOl For L = 30 we have 
r ~ 0.6 /iS, and for L = 40, r ~ 0.1 milliseconds. Thus, there is an angular momentum 
window that may yield mean lifes compatible with the experimental findings. Increasing L 
much further would produce too a distant Mg+Mg system which would correspond to a two 
independent Mg impurities in a drop. 

To check whether this foam-like structure of the Mg aggregate also appears in '^He drops, 
we have carried out calculations for (Mg+Mg) @^Heiooo using the same density functional 



as in Ref. 
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and the method presented in this Section. Figure [TT] shows the energy of the 
(Mg+Mg)@^Heiooo complex as a function of d. For distances smaller than some 8.6 A we 
have found that the system has a tendency to collapse into a dimer -physically unreachable 
from our starting point, Eq. (12 7p . We are led to conclude that there is no barrier in the case 
of liquid ^He. The configuration corresponding to the closest d we have calculated is shown 
in Fig. [H 

We are now in the position to determine the effect of these weakly-bound systems on 
the LIF and R2PI experiments on ^He droplets containing more than one Mg atom. Notice 
that the bottom right panel of Fig. [8] shows that the helium bubbles have a non-zero static 
quadrupole moment, whereas they are spherically symmetric for one single Mg atom in 
the drop. It is precisely the existence of this static quadrupole moment that causes an 
additional separation between the S and 11 spectral components in the absorption spectra, 
which results, as a consequence of the broadening of each line, in an double-peak structure of 
the computed spectra, in semi- quantitative agreement with the experimental data. Details 
of our calculation are given in Appendix C, see Eq. ( JC6I) . Figure [I3] shows such two- 
peak structure corresponding to the specularly symmetric configuration displayed in Fig. 
IHl and indicates that the 282 nm structure observed in the experiments may be attributed 
to the distortion produced by neighbour Mg bubbles; these bubbles contribute incoherently 
to the absorption spectrum, and the relative intensity of the 282 and 279 nm peaks might 
refiect the different population of drops doped with one and two Mg atoms, since those 
hosting a compact cluster (dimer, trimer, etc), would not yield an absorption signal in the 
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neighborhood of the monomer 3s3p ^Pi ^- 3s^ ^Sq transition. Notice that a large static 
distortion of the hehum bubble could also arise if the Mg atom were in a shallow dimple at 
the drop surface, but Fig. E] discards this possibility. 

We finally note that we have not considered in our work another source for an additional 
splitting of the spectral lines of a Mg atom in the field produced by a neighboring one, 
i.e. the resonant dipole-dipole interaction occurring during the electronic excitation. This 
effect could in principle lead to an additional (but probably small, compared with the effect 
discussed here) splitting of the calculated hues. An accurate determination of the dipole 
moment is required for a proper inclusion of this effect, which is beyond the scope of the 
present paper. 

VI. SUMMARY 

We have obtained, within DFT, the structure of ^He droplets doped with Mg atoms and 
have discussed in detail the magnesium solvation properties. In agreement with previous 
DMC calculations,-!^ we have found that Mg is not fully solvated in small ^He drops, whereas 
it becomes fully solvated in large droplet. 

As a consequence of its interaction with the helium environment, it turns out that magne- 
sium is radially quite delocalized inside the droplets. This large delocalization provides a way 
to reconcile two contradictory results on the solution of one Mg atom in a ^He drop, namely 
center localization (LIE and R2PI experiments^"^), and surface localization (electron- impact 
ionization experiments^^) . 

We have calculated the absorption spectrum of magnesium in the vicinity of the 3s3p 
^Pi ^- 3s^ ^So transition. For the large Mg@^Heiooo droplet, where Mg is fully solvated, 
we reproduce the more intense component of the absorption line found by LIE and R2PI 
experiments in large drops and in liquid helium. This agreement is only achieved when 
homegeneous broadening due to the coupling of the dipole excitation with the quadrupole 
deformations of the helium bubble are fully taken into account. This coupling is naturally 
included in Quantum Monte Carlo simulations of the absorption spectrum,— i^ whereby one 
takes advantage of the inherent fluctuations present in these simulations. These fluctuations 
are the full quantal equivalent of the dynamical distortions of the helium bubble we have 
introduced for the description of homogeneous broadening. An alternative method to include 
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shape fluctuations within DFT has been proposed and apphed to the case of Cs in bulk hquid 
hehum.— It would be interesting to adapt this method to the drop geometry, since it is not 
simple to handle dynamical bubble distortions in a non-spherical environment, or in ^He 
drops. 

To explain the origin of the low energy peak in the absorption line and conflrm the 
likely existence of soft, "foam" -like structure of Mg aggregates in ^He drops as proposed by 
Przystawik et al,- we have addressed the properties of two Mg atoms in ^Heiooo and have 
found that indeed, Mg atoms are kept apart by the presence of helium atoms that prevent 
the formation of a compact Mg cluster. We have estimated that the height of the energy 
barrier for the formation of the Mg dimer in ^He drops is ~ 2 — 3 K, which should be enough, 
at the droplet experimental temperature of 0.4 K, to guarantee a relatively long lifetime to 
these weakly-bound Mg aggregates. We predict that, contrarily, Mg atoms adsorbed in ^He 
droplets do not form such metastable states. 

The presence of neighboring Mg atoms in these structures induces a static quadrupole 
deformation in the helium bubble accomodating a Mg atom. As a consequence, the dipole 
absorption line around the 3s3p ^Pi ^- 3s^ ^Sq transition splits. We attribute to this static 
quadrupole moment the origin of the low energy peak in the absorption line, and conflrm 
the suggestion made by the Rostock group^ that the splitting of the absorption line, rather 
than being due to a dynamical (Jahn- Teller) deformation of the helium bubble, is due to 
the presence of more than one magnesium atom in the same droplet. 

Our previous study on Ca doped helium drops^ and the present work show that DFT 
is able to quantitatively address the dipole absorption of dopants in ^He drops within the 
diatomics-in-molecules approach, provided the impurity-helium pair potentials are accu- 
rately determined. However, we want to point out that, while we have a consistent scenario 
that explains the results of LIE and R2PI experiments, the understanding of the electron- 



impact ionization experiment reported in Ref. Il3| still requires further analysis. Indeed, 
since Mg atoms may be in the bulk of the drop or just beneath the drop surface, the exper- 
imental ion yield curve should reflect both possibilities, whereas apparently it does not (see 



Fig. 2 of Ref. Il3l ). One possible explanation may be that for electron-impact experiments, 
a A^ = 10 000 drop is still small, so that the impurity is always close enough to the drop 
surface to make the Penning ionization process to prevail on the direct formation of a He^ 
ion. 
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APPENDIX A 

In this Appendix we obtain the energy of the doped drop up to second order in 
the deformation parameters and the hydrodynamic mass of the helium bubble. To this 
end, the helium order parameter and Mg wave function are expressed as ^I/(r, t) = 



p{r,t)exp[i^^S{r,t)] and $(r,t) = |$(r, t)| exp[i^^v9(r, t)], respectively. Neglecting 
the velocity-dependent terms of the Orsay-Trento functional that mimic backflow effects,— 
the total energy of the system is written as 

E =^mHe J dh p{r,t)\VS{r,t)\'' + ^mMg J dhmr,t)\^\V^{r,t)\^ 

+ -^ rdh\V^/^{^)\' + -^ [dh MHv,t)\\'+ fdh£{p) 
2 rriHe J ^ 2 niMg J J 

+ J Jdhd'r'mr,t)\'VMg-He{\r-r'\)p{r',t) , (Al) 

where the functions p(|$|) and S{(p) fulfill the continuity equations 

-^p(r,t) = V[p(r,t)V^(r,t)] (A2) 

-||$(r,t)r = V[|$(r,t)rVv.(r,t)] (A3) 

that allow to identify S{r,t) and v^(r, t) as velocity field potentials, and the first two terms 
in Eq. (lAll) as a collective kinetic energy, whose density we denote as t[p, S, |$|, v^]. Thus, 

E = T + V = J dh{t[p,S,m,v]+v[p,\m . (A4) 
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In the adiabatic approximation, the dynamics of the system requires the following steps: 
i) introduce a set of collective variables (or deformation parameters) [a;(t)]) that define the 
helium density, p(r, t) = p(r, [Q;(t)]); ii) for each helium configuration defined by [a], solve the 
time-independent Schrodinger equation obeyed by |$(r, [a])\; iii) obtain the potential surface 
V[p, |<l>|] by computing the static energy for each configuration; iv) determine the velocity 
field potentials S'(r, [tt(t)]) and (/^(r, [«(t)]) by solving the continuity equations; v) compute 
the collective kinetic energy to obtain the hydrodynamic mass, and vi) solve the equation 
of motion associated with the effective Hamiltonian written as a function the deformation 
parameters. 

We aim to describe harmonic deformations of a spherical helium bubble created by an 
impurity in the ground state, and have to determine the helium density p(r, [«(t)]) resulting 
from a change in the radial distance to the center of the spherical bubble induced by the 
[ait)] parameters: 

oo A 

r^r + Y: E c^x,it)Z,,ir) , (A5) 

\=0 fi=-\ 

where the A = 1 deformation is now introduced to allow for displacements of the bubble. We 
recall that the real spherical harmonics have been normalized as (Zx^lZy^') = 2X+T^AA'^/i/i'; 
with ^oo(^) = 1) ZiQ{r) = cos9, etc. The breathing mode corresponds to A = 0, an 
infinitesimal translation to A = 1 (provided the fiuid is incompressible), and a quadrupolar 
deformation to A = 2. To first order, the density can be written as 

oo A 

p(r,t) =. Po(r) +p'o(r) ^ ^ a,^(t)Z,^(f) , (A6) 

A=0^=-A 

where from now on, the prime will denote the derivative of the function with respect to its 
argument. 

1. Impurity w^ave function 

To first order, the wave function |$(r, [a])\ is written as 

oo A 

|$(r,t)| ^ Mr) + E E «A,(t)<)(r) . (A7) 

A=0 ^=-A 
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The amplitudes <I>L (r) are determined in first order perturbation theory from the multipole 



expansion of the impurity-hehum pair potential 

ldh'p',ir')Z,,ir')Vx-Hei\r-r'\) 

J dh'p',{r')Z,,{r') J: \/i_^,(r, r')Z,,{r)Z, 



<^(r) 



x^r) 



= ui'Hr)ZUr 



Am I 



(A8) 



r(l), 



which defines U^ (r). We obtain 



«^?j(r) 



r(i). 



2^ $„£(r)Z^„(r) 



nim 



r(i)l 



„ ZA + i to ~ CnA^ 



^A^l 



(A9) 



that shows that actually, $U is /i independent. 



Once we have obtained the wave function, we can compute the energy surface K[p, |$|]. 
Since we describe deformations around a spherically symmetric ground state, the first order 
term vanishes, and the derivative g\ = d'^v/dax^day^i is also spherically symmetric. We 
can evaluate the second order contribution to the collective potential energy as 

1 r ., d^v 



\/(2) 



t^,,2J ' daxf^dax'^^' 



Xfi X'fi 

2-jr 



ax,,Zx^J,{r)ax'^J,' Zx' ^J.'{r) 



XlM 



2A + 1 7o 



-E 



27r 



po.vo 
drr^gx{r)al^ 

A 



2A + 1 



i^r^E 



a 



Xfi 1 



(AlO) 



A=0 ' M=-A 

which defines -E]^ . The parameters [a] are the dynamical variables that describe the evolu- 
tion of the system. 



2. Velocity field potentials 



Introducing the expansion S{r,t) = J2xti'^Xfi{t)Sx{r)Zx^{r), where the dot denotes the 
time-derivative, the continuity equation for the liquid helium is, to first order. 



-}_^ax, Zx, ^7 - E«A, Zx, | — ^ 



Po 



d^^A 2d^A A(A + 1)~, 



f-A/yt^ /y Cl^ 



s^ 



■ (All) 
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Hence, 



dpo dpo d^A 



dr 



dr dr 



+ Po 



d^S 



H ; n JX 



f-\ry'^ ly (H 'J^ ly ^ 



(A12) 



When r — * oo, the density vanishes for a drop, and approaches pb for the hquid. In the later 
case, Eq. (IA12p reduces to the radial part of the Laplace equation 







d^^A , 2d5A A(A + 1);, 



dr^ 



+ - 



whose general solution is 



Sx(r) 



r dr 



A, r^ + ""' 



Sx 



A+l 



(A13) 



(A14) 



59. Let r,- be the first 



We have solved Eq. (lA12p adapting the method proposed in Ref. 
point where po{r) is significantly different from zero [po{f) = for r < ri\. At this point, 
Eq. f lAT2|) implies that 



dpo 



dr 



dpo d^A 



dr dr 



which determines the boundary condition at r^ 

d^A 



dr 



-1 



(A15) 



(A16) 



For the liquid, the other boundary condition is that when r — >■ oo, the solution behaves as 
in Eq. ( 1A14I) with Aa = to have a physically acceptable solution. If the bubble has a sharp 
surface of radius rj and the liquid is uniform, 5'A(r) is completely determined by Eq. ( JA14I1 
and the velocity field potential that fulfills Eq. (IA16P corresponds to the coefficients 



Ax 
Bx 







,A+2 



A + l 



(A17) 



To find the velocity field potential in a large drop, we have defined a radial distance r^, 
far from the bubble and from the drop surface at -R1/2, around which on one may consider 
that the density is that of the liquid. Starting from r = r^ with the liquid solution fixed 
by the coefficients given in Eq. flA17p . we have integrated inwards Eq. flA12|) . finding the 
solutions S^^^{r) and 5'A(r) that correspond, respectively, to the non-homegeneous and to 
the homegeneous differential equation that results by setting to zero the left hand side of 
Eq. (lA12p . The general solution that satisfy the boundary condition Eq. (1A16P is obtained 



as 



S,{r) = ST{r) + CSl{r) 



(A18) 
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with 



l _l_ '^^ 



dr 



c 



dlh 



(A19) 



ri 



The field ly? is analogously obtained after introducing the expansion (/?(r, t) = 
J2xfj.<^Xfj.{'t)^\{f^)^\fi{'^)- We have assumed that the wave function of the impurity in the 
ground state is a Gaussian $o(^) = ^exp {—f3r^) whose shape has been determined by fit- 
ting it to the actual wave function, and have introduce a cutoff distance Vg such that safely 
$o(^) = ii r > Tg. We then obtain the following differential equation to determine ipx: 



A/3r = -Aj3r 



d^x , d^A , 2d^A A(A+1). 



+ 



dr ' dr^ 

The appropriate boundary conditions are 

d(^A 



+ 



r dr 



Vx ■ 



(A20) 



dr 

''0 

^a(O) 







(A21) 



and the solution is analytical but very involved (we have found it by using the Mathematica 
packpage). The dipole mode is the only exception; if A = 1 we have ipi = — r. 

3. Kinetic energy 

Once the velocity fields S{r,t) and v^(r,t) have been determined, the collective kinetic 
energy can be easily calculated to second order in the collective parameters: 

T = Jdht[p,S,\<^\,v] 

^niHe J dh p{r,t)\VS{r,t)\^ + ^mMg J dhmr,t)\'\Vv{T,t)\' 
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ex 
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^ An 



A=0 



2A + 1 2 



rriHe / dr r^ po{r) 



'dSA ^ A(A + 1) ^, 
dr / r^ 



^'^Mg J dr r^ |$o(^)P ( 



dr / r^ 



E« 
fj,=—\ 



A/1 



1 oo A 

- ^E^A* E «l 

^ A=0 fj.=-X 

which defines the hydrodynamic mass M^ for each A mode. 



(A22) 
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Using the continuity equations and the Gauss theorem, one can find an alternative ex- 
pression for M^ 

4:71 



Ml 



^y^ [rriHe j dr r^ p'^{r) Sx{r) + mug j dr r^ 2%{r)%{r)(px{r)^ . (A23) 



We have checked that both expressions yield the same values for M^^ which constitutes a 
test on the numerical accuracy of the method. It is easy to see from Eq. flA22p that in 



59. 



bulk liquid helium, the A = 1 hydrodynamic mass coincides with that given in Ref. 
Using that ipi = — r, it is also easy to check from the above expressions that the impurity 
contribution to the A = 1 hydrodynamic mass is just the bare mass of the Mg atom. 

The sum of Eqs. (JAIOI) and flA22p represents the Hamiltonian of a set of uncoupled 
harmonic oscillators whose frequency only depends on A. 

APPENDIX B 

In this Appendix we work out in detail the expressions we have used to describe the 
homogeneous broadening of the absorption line. We consider a spherical ground state defined 
by a helium density po{r) and impurity wave function ^o{r), both modified by the action 
of the breathing mode defined in Eq. ( 12T1) . namely p(r, ao) = po{r + oq) AI^~\ and ^(r, ao)- 
The computation of the spectra for a given ao can be carried out starting from Eq. ( ITTl) 



I{u,ao)=A7rY,Jdr \r <l>{r,ao)\'6[u + u{ao) - Vr{r,ao)/h] , (Bl) 

i 

where hu^a) = e{ao) is the impurity eigenenergy and V^^{r,ao) are the PES defined by 
Eq. (IT3ll . where [a] reduces to ao- Notice that ao is not introduced perturbatively; it is 
unnecessary since this mode does not break the spherical symmetry. Next, we perturbatively 
introduce the modes with A > 2 to first order; Eq. fIBll) becomes 

I{uj,[a]) ~ I{uj,ao) 

oo A 



-47r J dr \r <l>(r, ao)|'<5'[^ + ^(ao) - Vr{r, ao)/h]^e\^{r)^ ax^ , (B2) 



where e\^(r) = dVf^{r)/daxfj,\ia^ =o- "^^^ ^^^^ integral is zero due to the orthogonality of 
the spherical harmonics. To evalue the second integral, we expand e\^(r) as a power series 
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of r. Taken into account that the PES have a stationary point at r = due to the spherical 
geometry -the first order term is zero- we can safely stop the expansion at the zeroth order 
term, since the wave function $(r, a^) is very narrow 



I{uj,[a]) c^I{uj,ao) 

oo A 



-4^EE E fdr\rHr,ao)mu~u;{ao)-Vr{r,ao)/h](-e\^{0)a,,) . (B3) 

i A=2m=-A ^ ^ 

This equation may be interpreted as the expansion to first order of a shift in uj, so it can be 
written as 

oo A 



I{u,[a]) ~ 47r^ / dr |r $(r, ao)r<^ 

i 



uj + cj(ao) 



h 



yr(r,«o) + E E ev(o)«A, 

A=2^=-A 



oo A 

^+E E ev(0)aAMA,«o 



(B4) 



A=2/x=-A 

Thus the eigenvalues e\^(0) are related to the diagonalization of the expansion of Uij{r, [a]) 
defined in Eq. flT2l) . Writing this matrix equation as a function of the real spherical har- 
monics 



f/(r,[«])=y'dVp(r' + r,H) 



/ 



-[Vi(r')+2Vn(r')]^o(r) 



1 
1 
1 



(B5) 



+ ^[V^{r')-Vn{r')] 



/-^Z2o(r)+Z22(r) Z2_2(r) Z^iif) M 

Z2-2{f) -^Z2o(r) - Z22{r) Z2-i{r) 



\ 



^21 (r) 



Z2-i{f-) ^^2o(r) I J 



expanding p(r' + r, [a]) to first order in [q;]a>2, and evaluating the first order contribution 
at r = we obtain 

t/(r,H) ^ f/(r,ao) + E E 

A=2/x=-A 



2A + 1 7 ^'^'^'^P'('^''«o)-^[Vs(r') - Vu{r')] 



I — l^S2x6of, + 62x62,, 



X 



6o\6-': 



6o\6^ \ 



V 



"Am 



(B6) 



^2A(5-2m --^52A50At - 62\62^_i 62x6-!^ 

62x6ipi 62x6-1^, -^62x601, 1 

This shows that, to first order, only quadrupolar deformations are coupled to the dipole 
electronic transition, and that its effect is a shift of the spectral line, as shown in Eq. (1B4I) . 
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At variance with the approach of Ref. |29|, where the above matrix is approximated by its 
diagonal expression, implying that only the fi = and 2 components of the quadrupole 
deformation are considered, our approach incorporates all five components. 

The relation between the eigenvalues Aj(r, [a]) of U{r, [a]) and the coefficients e^ (0) is 

2 






Finally, the total spectrum is written as 



I{uj) oc An daod^ a2\ip{ao, 02)1"^ I {uj, [a]) 

~ 47r / daod^a2\^{ao,a2)\'^Yl / '^^ 1^ ^i''",ao) 
x6 u + t^f (ao) - T^ii^^ N) 



(B7) 



(B8) 



APPENDIX C 



In this Appendix we consider that the doped drop is cylindrically symmetric. We expand 
the helium density and ground state wave function into spherical harmonics with /x = 



A=0 

00 

<l>(r) = E7A(r)ZAo(r), 



(CI) 



A=0 



where /3a = ^ J dhp{r)Zxo{r) and px{r) = ^ J dnp{r)Zxo{f), with analogous defini- 
tions for $A and 7a. If Po ^ Px>o (^ 7o ^ 7a>o)7 we can compute the line shape to first 
order in /5a>o (7a>o); in analogy with Eq. flB2p we write 



-4^E E y dr \r 7o$o(r)p5'[a; - ^(ao) - Vr{r, Po)mj-eMPx 



(C2) 



i A=l' 



with e\o defined now by the potential matrix 



4.7]- r 

Uir,[(3]) ^ f/(r,/5o) + — ydr'r'2p2(r')[rs(r')-Vn(r')] 



-1 


o\ 





-1 





2) 



P2 ■ (C3) 
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Introducing the shape deformations defined in Eq. (lASp . to first order this matrix becomes 



[/(r, [/?],[«]) ^ f/(r,/3o,«o) + 



5^3 



dr' r' Po p'^ir', a^) [Vs(r') - Vnir')] 



I -^«20 + "22 a2-2 



X 



^21 



\ 



V 



^2-2 
"21 



^20 — ^22 <^2-l 



v^ 



(C4) 



"2-1 ^'^20 / 



where we have defined 020 = "20 + P2C , with 

/dr'r'2p2(r')[Vs(r')-Vh(r')] 

C = ^ . (C5) 

jdr'r''PoP'o{r\ao)[V^{r')-Vn{r')] 

These equations show that the computation of the fine shape for this geometry is as in the 
spherical case but with a shift in the 020 parameter. 
The dipole absorption spectrum is finally obtained as 

|2 



I{ijj) (X An dao d^a2|'?/'(ao5 "2)1^ X! / '^''" I'" '^('^'Q^o) 



x6 



u; + u;f (ao) - T^iir, [a]) 



(C6) 



with a'l = (^20 - P2CY + ^22 + "2-2 + "21 + "Li- 
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TABLE I: Atomic shift Acu of Mg@^Heiooo (-^1/2 = 22.2 A) as a function of the average distance 
between the magnesium atom and the center of mass of the ^Heiooo nioiety. Also indicated is the 
corresponding wavelength A. The value of the transition energy in the gas phase is 35 051 cm~^j^ 
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FIG. 1: (Color online 
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X^T, Mg-He pair potentials used in this work: 
19l ; circles connected with a dotted line, from Ref. 
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FIG. 2: Energy of the Mg atom as a function of the number of atoms in the drop, obtained using 
the Mg-He potential of Ref. 19 (squares). The values given in Ref. [tI are also displayed (dots). 
The lines are drawn to guide the eye. 
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FIG. 3: (Color online) Bottom panel: total energy (K) of Mg®*^! 



n 



e^o as a function of Zq (A) 



2d (circles). The energies are 



obtained using the Mg-He potential of Ref. Il9| (squares) and of Ref. 
referred to their equilibrium values, -157.0 K and -153.8 K, respectively. The vertical line locates 
the drop surface at i?i/2 = tqA^^'^, with rg = 2.22 A. The horizontal line has been drawn 0.4 K 
above the equilibrium energy. Top panel: probability densities for the configurations displayed in 
the bottom panel; the single peak distribution corresponds to the Mg-He interaction of Ref. l20. 
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FIG. 4: (Color online) Helium densit y p lots of the Mg@ Heso droplet in the y = plane obtained 



using the Mg-He potential of Ref. |19|. From top to bottom and left to right, the Zq values 
correspond to the stationary points displayed in Fig. [3l namely 0, 2, 4, and 6 A, respectively. The 
brighter regions are the higher density ones. 
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FIG. 5: (Color online) Helium density profiles of the Mg@^He5o droplet along the z axis obtained 
using the Mg-He potential of Ref. Il9|. The Zq value is indicated in each panel. Dashed lines, pure 
drops; solid lines, doped drops. 
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FIG. 6: (Color online) Bottom panel: total energy (K) of Mg@^Hein nn a s a function of Zq (A) 
obtained using the Mg-He potential of Ref. 19 (squares) and of Ref. 20l (circles). The energies 
are referred to their equilibrium values, -5482.0 K and -5476.8 K, respectively. The vertical lines 
roughly delimit the drop surface region, conventionally defined as the radial distance between the 
points where the density equals O.l/jfe and 0.9pfe, being p\, = 0.0218 A^^ the bulk liquid density. The 
horizontal line has been drawn 0.4 K above the equilibrium energy. Top panel: total energy (K) of 
Mg@^HeAr with N = 1000 and 2000 (vertical right scale) as a function of Zq (A) obtained using 



the Mg-He potential of Ref. Il9|. The energies are referred to their equilibrium values, —5482.0 K 
and —11629.8 K, respectively, and the distances (horizontal scale) are referred to the R1/2 radius; 
also shown are the corresponding probability densities (vertical left scale): dot (dot-dash) line, 
N = 1000(2000). The solid line represents the probability density of the A^ = 10 000 drop. 
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FIG. 7: (Color online) Total absorption spectrum of one Mg atom attached to ^Heiooo in the vicinity 
of the 3s3p ^Pi <— 3s^ ^Sq transition. The line has been decomposed into its two 11 and one E 
components, the former one is the higher frequency transition. The starred vertical line represents 
the gas-phase transition. The experimental curve, adapted for Ref. 8, has been vertically offset 
for clarity. Also shown is the absorption spectrum obtained neglecting homogeneous broadening 
(hatched region) 



44 



20 



N 



-20 




20 



-20 




20 



°S 



-20 




-20 

x(A) 



20 




-20 20 

x(A) 



FIG. 8: (Color online) From top to bottom and left to right, (Mg+Mg)@ Heiooo metastable con- 
figurations for Mg-Mg interatomic distances d = 18.5 A, 12.9 A, and 9.3 A, and total energies 
-5567.8 K, -5573.9 K, and -5580.3 K, respectively. The bottom right panel shows the specularly 
symmetric configuration at d = 9.5 A with total energy —5581.4 K. The brigther regions are the 
higher density ones. 
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FIG. 9: (Color online) Helium density profiles of the (Mg+Mg)@ Heiooo symmetric configuration 
{d = 9.5 A) along the z-axis (solid line) and the x- or y-axis (dashed line). 
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FIG. 10: (Color online) Energy (K) of the Mg+Mg@^Heiooo system as a function of the Mg-Mg 
distance (A). The energies have been referred to that of the metastable equilibrium configuration 
(local minimum) at L = 0. The lines have been obtained by a cubic spline of the actual calculations. 
From bottom to top, the curves correspond to L = to 50h in Wh steps. 
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FIG. 11: (Color online) Energy (K) of the Mg+Mg@'^Heiooo system as a function of the Mg-Mg 
distance (A). The Hne has been drawn to guide the eye. 
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FIG. 12: (Color online) Top panel: helium density profiles of the (Mg+Mg)@^Heiooo complex at 
d = 8.6 A along the z-axis (solid line) and the x- or y-axis (dashed line). Bottom panel: equidensity 
lines corresponding to the same configuration. The brigther regions are the higher density ones. 
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FIG. 13: (Color online) Total absorption spectrum of one Mg atom attached to ^Heiooo in the 
vicinity of the 3s3p ^Pi <— 3,s^ ^Sq transition in the distorted environment created by the presence 
of another Mg atom. The line has been decomposed into its two 11 and one S components, the 
former one is the higher frequency transition. The starred vertical line represents the gas-phase 
transition. The experimental curve, adapted for Ref. 8, has been vertically offset for clarity. Also 
shown is the absorption spectrum obtained neglecting homogeneous broadening (hatched region) 
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